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ABSTRACT 



Aims. We investigate the distribution of the luminous and the dark matter components in the isolated ellipticals NGC 7052 and 

NGC 7785, which are embedded in an emitting hot gas halo, by means of relevant X-ray and photometric data. 

Methods. To calculate the dark matter distribution in these rare objects, we performed an improved X-ray analysis of the XMM- 

Newton data of NGC 7785, and we used former results based on Chandra data of NGC 7052. For each object we also derived the 

stellar spheroid length scale from the surface photometry and the spheroid stellar mass from an analysis of the galaxy spectral energy 

distribution. 

Results. We find that a dark matter component is present in these objects. It is subdominant and mixed with the luminous matter 
inside the optical region half-light radius wide, while it dominates the gravitational potential at outer radii. On the whole, the dark 
halo structure is very similar to that found around spirals of comparable luminosity and it is well reproduced by a Burkert halo, while 
a Sersic spheroid accounts well for the baryonic component. 

Key words. Galaxies: elliptical and lenticular, cD - X-rays: galaxies - Galaxies: photometry - Dark Matter 



1. Introduction 

The dynamics of galaxies is influenced, and often dominated, 
by non-radiating matter, which reveals itself only through a 
gravitational interaction with the luminous matter. From the 
study of a large number of disk systems, a one-to-one rela- 
tion between the (luminous) central object and the massive dark 
halo that envelopes it h as emerged (see IShankar et al.l l2006t 
iMandelbaum et all l2Q09h . The paradigm is that every galaxy 
lies within a self-gravitating dark halo of size R w i r and mass 
M v j r : more specifically, disks of sizes R opt and luminous masses 
M+ are embedded in dark halos of virial radii /? v j r ~ 15/? opt 

and masses M vir = 3 x 10 12 ( ~^\ .. )° 4 Af B , where M Q is 
the solar mass (jPersic et aflll996t iKormendv & Freemanll2004t 
ISalucci et al.ll2007t lOh et al.ll2008h iDonato et alj|2009l 

D o ellipticals follow a s i milar scenario (see iBertola et alj 

I1993L ISalucci & Persidfl997t iTortora et al.ll2009h ? The answer 
would be of crucial importance for understanding how galaxies 
form and perhaps even for the nature of the dark matter (DM) 
itself. Unfortunately, evidence on the properties of the mass dis- 
tribution in these objects is presently not conclusive. 

Indeed, to derive the mass distribution from stellar kine- 
matics, as is done in spirals where rotation balances gravity, is 
very challenging. Ellipticals are pressure-supported stellar sys- 
tems whose orbital structure may also involve some angular mo- 
mentum content, a degree of triaxiality, and velocity dispersion 
anisotropies. As is well known, the derivation of the mass distri- 
bution from the stellar motions in these systems through Jeans 
equation and its higher moments is quite complex. It is difficult 
to obtain the galaxy gravitational potential unambiguously and, 
following this, the mass distribution. Furthermore, in these sys- 
tems the luminous matter is very concentrated toward the cen- 
ter, where it dominates the gravitational potential. Therefore, 



the motions of the main baryonic component in ellipticals are 
generally a limited tracer of the distribution of the dark matter, 
which mostly lies outside the stellar spheroid (e.g., Kronawitter 
et al. 2000). Only recent dynamical measurements reached the 
DM dominated regions (e.g., Thomas et al. 2009 and references 
therein; Pu et al. 2010). 

Alternative mass tracers (planetary nebulae, ionized globular 
clusters, and neutral disks) from baryonic components of neg- 
ligible mass a re often a better pro be of the total gravitational 
potential (see lRomanowskvll2010l and references therein). In 
particular, the X-ray emitting hot gas is known to extend out 
to very large radii and to trace the gravitat ional potential. This 
method yields accurat e mass profiles (e.g . , iForman etaj] ll98a 
IP' Sullivan et al. 200^ IZhang et al.lr2007l) and it is very robust 
because it provides reliable mass profiles also for highly dis- 
turbed objects where the hydrostatic equilibrium is far from be- 
ing established dJohnson et al. 2009|). 

This method works best for isolated systems and has the rare 
advantage that it can probe the gravitational potential from al- 
most the galaxy center out to the halo virial radius. Moreover, to 
obtain mass profiles of objects that never resided in a cluster or 
group resolves the question of whether ellipticals possess their 
own DM halos, or "build them" from the cluster DM, and it al- 
lows us to investigate the environmental dependence of the DM 
distribution in galaxies. 

2. X-ray emitting halo ellipticals 

A regular X-ray morphology is commonly assumed to indicate 
approximate hydrostatic equilibrium, and the absence of nearby 
companions leaves the X-ray halo unperturbed. The gravitating 
mass inside a radius r, M(r), can be derived from the X-ray flux 
under the assumption that the emitting gas is in hydrostatic equi- 
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librium. As is well known, from the gas den sity and temperature 
profiles (e.g jFabricant & Goresteiril 19831) we obtain: 



„„, s kT g (r)r(dlog p g d\ogT g (r)\ 

M{< r) = — — + — 

G/i»jp V a log r a log r ) 



(1) 



where T„ is the gas temperature at radius r, p g is the gas density, 
k is the Boltzmann constant, G is the gravitational constant, fi is 
the mean molecular we ight that we set at the value of p = 0.62 
dEttori & Fabian ||2006). and m p is the mass of the proton. 

The observed X-ray surface brightness, Z g is usually fitted 
by ayS-model (Cavaliere & Fusco-Femiano 1976): 



E g oc [1 + (r/r c ) 2 ]-^ 2 
that yields 

p g (x[l+(r/r c ) 2 ]- 3 ^ 2 , 



(2) 



(3) 



with the core radius r c , and the slope B as free parameters. For 
an isothermal gas distribution: dlogT g (r)/dlogr = 0, therefore, 
we have, in a spherical case, from the above equations: 

Wr c ) 3 



M(r) 



[l+W'-c) 2 ]' 



(4) 



An isothermal distribution is often found in ellipticals; how- 
ever, there are cases in which the temperature has a profile that 
must be considered in Eq. (1), although its contribution is al- 
most al ways much smaller than that from the gas density gradi- 
ent (e.g. lFukazawa et aT]|2006l) . In the present work, we have a 
very low spatial resolution temperature measurement, therefore 
we are forced to assume that the temperature in the region under 
analysis, varies linearly with radius. An assumption that, how- 
ever, agrees well with most of the high-resolution temperature 
profiles. Then, Eq. (4) or, when the temperature gradient is dif- 
ferent from zero, Eq. (1) yields the gravitating mass as a function 
of the radius in the region where the X-ray halo emits. 



3. The isolated ellipticals NGC 7052 and NGC 7785 

The two elliptical galaxies we study here are members of a 
sample of 43 bright isolated galaxies dFocardi & Kelm 201 ll) 
select ed from the Updated Zwicky Catalog ( UZC, Falco et al. 
by means of an adapted version of the iFocardi & Kelm 



neighbor-search code. Seven of those are early-type ob- 
jects, for four of them results of the X-ray anal ysis of propri- 
etary and archival data have been published by iMemola et alJ 
d2009l) . and finally, for two of them (NGC 7052 and NGC 7785) 
the available X-ray statistics allows us to derive the dark matter 
content and distribution. 

We recall the criteria to select the sample of truly isolated 
galaxies: a) absolute B luminosity higher than 1.3x 10 l0 /i^ Lb©, 
where Lb g is the fi-band solar luminosity; b) recessional velocity 
in the range [2500-5000] km s" 1 ; c) galactic latitude (fo n| > 15°); 
and d) no companion galaxies down to 2 magnitudes fainter than 
the target galaxy within a circular radius of l.3hj} Mpc and at 
velocity distance of 1000 km s -1 . These ellipticals are selected in 
an objective way and their isolation is granted at a typical clus- 
ter/group scale. Moreover, NG C 7785 is present in the AMIGA 
database of isolated galaxies ( Verdes-Montenegro et al.l 120051) 
and is also part of the lSmith et al .1 (12004 sample. 

The importance of selecting objects in total isolation was 
discussed in the introduction: dark matter halos of cluster el- 
lipticals merge with the cluster halo at a scale of ~ 300 kpc, a 



Table 1. Physical properties of the galaxies. Columns: (1) galaxy 
name; (2) core radius of the best-fit /?-model; (3) slope B; (4) 
absolute B-band magnitude; (5) B-band magnitude; (6) A^-band 
magnitude, 2MASS total magnitude; (7) scale in kpc/", which 
also provides the distance assumed. B- and K s -bwd magnitudes 
are from LEDA or NED. 



Source 


r c 


B 


M B 


rriB TC 


K 


kpc/" 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 7052 


1.11" 


0.48 


-21.46 


12.73 


8.57 


0.323 


NGC 7785 


35" 


0.95 


-21.38 


12.34 


8.45 


0.247 



scale comparable with that of their virial radii. This complicates 
the dynamical analysis, and more importantly, it implies that the 
present-day halos around cluster galaxies are not the cosmolog- 
ical dark matter halos anymore, indeed, ellipticals at the center 
of clusters are f ound to have more DM tha n ellipticals in low- 
density regions dNagino & Matsus hita 2009). The very low den- 
sity of the environment of isolated ellipticals makes them pre- 
cious laboratories for investigating the cosmological properties 
of DM halos around early-type galaxies. 

Note that NGC 7052 and NGC 7785, with L x < 2 x 10 41 
ergs -1 , and M TO t < 2 x 10 12 M Q (see Table 3 in Memola et 
al. 2009) do not appear to be fossil groups (e.g., Jones at al. 2003; 
Pompei et al. 2007) because they do not show their physical 
properties. Fossil groups (e.g., NGC 1 132, Mulchaey&Zabludoff 
1999) have, indeed, so far been observed with larger masses 
10° - 1O 14 M , and X-ray luminosity > 10 42 ergs -1 (see Fig. 8 
in Memola et al. 2009). 

4. X-ray observations 

4.1. X-ray data 

NGC 7052 and NGC 7785 have extended X-ray halos out to 
16 and 32 kpc, respectively, with hot gaseous mass of M aas = 
2.2 x 10 9 M o , and M gas = 4.6 x 1O 9 M . Chandra observations 
of NGC 7052 and XMM-Newton obse rvations of NGC 778 5 
have been presented and discussed in IMemola et al.l d2009l) . 
NGC 7052 was observed by the Chandra ACIS-S instrument 
in September 2002 for about 10 ks, and NGC 7785 was ob- 
served by the XMM-Newton EPIC detectors in June 2004 for 
about 15 ks (Prime Full Window mode a nd thin filter). For th e 
spectral analysis of NGC 7052 we refer to IMemola et all d2009l) . 

The net radial surface brightness profiles, centered at the X- 
ray peak of the EPIC-pn data in the 0.5-2.0 keV energy band, are 
well reproduced by a ff-mod el whose best-fit valu es of r c and B 
are given in TableQ](see also IMemola et al.ll2009l) . 

4.2. The X-ray spectral analysis 

We present another X-ray spectral analysis of the NGC 7785 
data. Thanks to the quality of the X-ray data and in view of the 
structure of Eq. (Q]), we aim to deeply investigate the tempera- 
ture profile of the X-ray gas emitting component of this isolated 
elliptical galaxy. 

Source counts were extracted for the XMM-Newton EPIC- 
MOSl,-MOS2, and -pn from two separate, but contiguous, re- 
gions: a 65"radius circle centered on the source, and a concentric 
annulus with 65"and 130"radii. The background spectrum was 
extracted from a circular source-free region of 90", at 3.8 arcmin 
distance from the target. The X-ray spectra were analyzed with 
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NGC 7785 NGC 7785 




Energy (keV) 01 0.5 1 2 5 

Energy (keV) 



Fig. 1. XMM-Newton EPIC-MOS1 (red), -MOS2 (green), and -pn (black) spectrum and x 2 behavior for NGC 7785. Left panel: 
extraction radius of 65"; Right panel: extraction annulus with 65"and 130"radii. The best-fit model (solid line) is a power-law plus 
a thermal component (see color figure online). 

Table 2. NGC 7785: summary of the XMM-Newton spectral results. The EPIC-pn net counts are derived from the regions (a circle 
and an annulus) of radii given in the second column. 
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cts 
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keV 


r 

fixed 


X 2 ldof 


Lmekal(0.5-2.0) 

erg s"" 1 


Lp O (2.0-10) 

erg s -1 


5.24 xlO 20 


65 


2687+59 
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307/314 


6.3 xlO 40 


6.9 xlO 39 


5.24 xlO 20 
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615±64 
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39/39 


9.4 xlO 39 


5.4 xlO 39 




Fig. 2. Two-dimensional 2MASS /<f s -band image decomposition of NGC 7785. From left to right: the original image, the model 
image, the residual image, and the surface brightness profile plot. The image surface brightness profile is marked by filled circles 
and the model with a solid line. The model-data residuals are also shown. 



the XSPEC package (version 12.6.0: lArnaudlll996l) . The source 
counts were stacked into energy bins such that each bin has a 
significance of at least 2cr, after background subtraction. The 
quoted errors on the best-fit parameters correspond to the 90% 
confidence level for one interesting parameter (i. e., Ax 2 = 2.71; 
lAvnil[T976T) . 



In both cases the best-fit function of the source counts to 
fit the hot gas includes a power-law with index fixed at T =2, 
which models the central X-ray binaries or AGN component, 
plus a thermal (mekal model) component with abundances fixed 
at 50% the solar values. For both components the appropriate 
Galactic hydrogen column density (Na) along the line of sight 
dDickev & Lockmanl Il990l) has been taken into account. We 
found (see Fig. Q] and Table H kT g = 0.32+Jj ^ keV in the outer 



ann ulus and kT s = 0.57 + » keV in the inner circle, in agreement 
with lMemola et alJ d2009f) . Hereafter, we assume r me k a i = T g . 

A best-fit for the outer region (x 2 /d.o.f. = 39/39) with two 
free parameters gives the kT = 0.32^2!!? keV value, which is in- 
compatible with kT = 0.56 + 0.02 keV that refers to the hypote- 
sis of constant temperature. Our result is fairly robust: a constant 
value of kT g out to 130 arcsec can be excluded at > 4<x level, and 
consequently a temperature gradient emerges that we fit linearly. 

In detail, the data are well reproduced by 

kT g (r) = [0.6 - 1.2 x 10~ 2 (r - 3)]keV, (5) 

with r in kpc. This result is not surpri sing and agrees with 
the gen eral behavior of isolated ellipticals: iNagino & Matsushita! 
(120091) find that these objects tend to show a negative tempera- 
ture gradient. Notice that this variation implies in Eq. (Q]i a mild 
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contribution from the temperature gradient term that it is smaller 
than the one from the gas density gradient. Therefore, the effect 
of the uncertainty in the temperature profile on the mass distri- 
bution of Eq. © is negligible. 

In NGC 7052 we cannot statistically support any evidence 
of gas temperatu re variation. F o r this galaxy, then, dT/dr = 
as claimed in iMemola et all d2009t) . A constant tempera- 
ture has also been observed in so me other isolated ellip ticals, 
e.g. NGC 57, NGC 7796, IC 1531 dO' Sullivan et alJl2007h . 



5. Structural parameters of the stellar spheroids 

NGC 7785 is known to display boxy isophotes dLauerl 19851) . and 
taking this into account, we re-derive its structural parameters 
from the /<f s -band image (see Fig.|2| taken from the Two Micron 
All Sky Survey (2MASS) atlas. The pixel size is 1 arcsec, and 
the point-spread function full-width half maximum is 2.9 arcsec. 
We used the two- dimensional galaxy fitting algorithm GALFIT 
dPeng et al.ll2002l) for fitting one Sersic component, which is al- 
lowed to have disky or boxy isophotes. The sky value is fixed to 
the value provided in the 2MASS image header. The half-light 
radius is found to be r e = 22.70" and the Sersic index « = 3.8, 
very close to the pre viously found de Vaucouleurs fit r e = 19.5" 
dBender et al.lll988h . The GALFIT boxiness parameter for the 
best-fit model is Co = 0.19, and we verified that the addition of 
this parameter in the case of this galaxy does not change the half- 
light radius or the Sersic index. If, instead of fixing the sky value, 
we let it be a free parameter, which is determined by GALFIT, 
both the r e and n values are 10% lower. This effectively sets the 
uncertainty of the fit, which is negligible for the results in the 
succeeding sections of this work. The performance of this fit can 
be seen from the residual image in Fig.|2](c). 

We also compared the surface brightness profile from the 
above fit with the surface brightness profile of the galaxy image 
constructed by means of the IRAF task ELLIPSE. The result is 
shown in Fig. [2] 

In NGC 7052 the photometry is simple and well studied 
bv lBenderetai]dl988l) : the distribution of the stars is well de- 
scribed by the usual de Vaucouleurs spheroid with half light 
radius of r e = 24". We confirmed this value with a two- 
dimensional image decomposition of a 2MASS ,ff s -band image 
performed as above. However, we omit the plots for brevity rea- 



6. Stellar spheroid mass in NGC 7052 and NGC 7785 

It is crucial to obtain an independent estimate of the mass of 
the stellar spheroid for these objects. We computed them by fit- 
ting the observed spectral energy distribution (SED) with syn- 
thetic spectra. We proceeded as follows: photometric data in 
B,V, and 2MASS J,H,K S bands were adopted from NED (ho- 
moge nized NED va l ues) a nd corrected for Galactic extinction 
from ISchlegel et al.l d!998l) . We used the total magnitudes in 
2MASS bands, which are obtained by integrating the radial 
surface brightness profile out to r tot = 87.1" and 82.2" for 
NGC 7052 and NGC 7785, respectively. Given that half-light 
radii for these two galaxies are r e = 24" and 22.7", integration 
out to nearly four times thes e radii provides a good estimate of 
the total galaxy magnitude dJarrett et al. I l2000h . The error esti- 
mate for the total magnitude is 0.02 mag for both galaxies. The 
uncertainty introduced by the extinction correction and possible 
systematic differences between bands is accounted for by con- 
sidering an additional 0.05 mag error in quadrature. Photometric 



Fig. 3. NGC 7052. Left: Model fit to photometry (see color fig- 
ure online). The best-fit spectral template is marked with points, 
the corresponding model magnitudes in the observed bands with 
circles. The observed magnitudes are indicated with filled red 
circles and error bars. The (black) bars in the wavelength direc- 
tion denote the bandwidth of respective filters. Right: The likeli- 
hood of the stellar mass estimates. 



- / 














a". 

v.* 



5000 10 4 1.5x10* SxlO 4 2.5xl0 4 10.! 

\ / K 

Fig. 4. Same as Fig.[3]for NGC 7785. 



data for NGC 7785 were available in the above five bands, while 
for NGC 7052 the V-band data were unavailable. 

The stellar mass is obtained from the mass-to-light ratio of 
the best-fitting model and the K s values (see Table 1). We fit 
the synthetic spectra to the observed spectral energy distribu- 
tion by using the modified version hyperzmass of the HyperZ 
code dBolzonella et alJ2000h . kindly provided by M. Bolzonella. 
The redshifts are fixed to z — 0.0156 and z = 0.0127, for 
NGC 7052 and NGC 7785, respectively. The adopted set of 
synthetic spectra consists of composite stellar populations with 
two bursts of star formation, b ased on simple stellar popula- 
tion models of Maraston] ( 2005 ) and constructed using the code 
GALAXEV dBruzual & Charlotll2003l) . The Salpeter initial mass 
function and solar metallicity are adopted. The strength of the 
more recent burst is 5% of that of the older one, and they are 
separated by a period varying in the range 0-11 Gyr. The adopted 
star-formation rate is consistent with what is typically found for 
elliptical galaxies and includes an old major burst at approxi- 
mately the halo-formation redshift and the presence of a weaker 
burst in the last several gigayears. 

For NGC 7052 the best-fit model (Fig. [3] filled black cir- 
cles) reproduces the data well (model magnitudes: open circles; 
data: filled red circles and error bars). It consists of two bursts of 
star formation, one 1 1.3 Gyr ago, and the second one, of relative 
strength 5%, still occurring. The resulting stellar spheroid spec- 

at 2<t fitting 



trophotometric mass is log(M'^ /i /M ) = 11.6^ ', 
uncertainty, then log(M'^ /i /Ls) 



For NGC 7785 the best-fit model reproduces the data well, 
in this case it points to a main burst 9.0 Gyr ago, followed by 
the second weak one 8.0 Gyr ago. For this object, a single burst 
of star formation about 9.0 Gyr ago would also be a satisfactory 
model. The stellar spheroid spectrophotometric mass is found 
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at 2cr fitting uncertainty, then 



to be \og(M' iph /M B ) 
log(M* ph /L B ) = 0.20!°;3. 

Figures [3] and [4] show the spectral template fitting (left) and 
the mass determination (right) for NGC 7052 and NGC 7785, re- 
spectively. In both cases a large number of different SFRs yield 
a stellar mass similar to the one we assume. Models in the like- 
lihood plot differ by the age of the two components and the red- 
dening, which is for simplicity assumed to be the same for both 
components. The mass estimates obtained by this method are 
very robust: they primarily depend on the normalization of the 
SED in the near-infrared bands, and are relatively insensitive t o 
the details of the star-formation history dMaraston et al.ll2010l) . 
In particular, changing the star-formation history to a single star- 
formation episode of varying length and functional form results 
in a maximum 10% change in the stellar mass (for the best-fit 
model) for NGC 7785 and 20% for NGC 7052; varying metal- 
licity in the range between half-solar and two solar metallicity 
leaves the mass unchanged for NGC 7785 and 30% smaller for 
NGC 7052; changing the relative strength of the second burst 
from 5% to 10% gives the same mass for NGC 7785 and a 10% 
lower mass for NGC 7052. Again, these modifications mainly in- 
fluence the implied age for the population, with a small effect on 
the deduced stellar mass, which is well within the quoted error 
bars. Using the Chabrier or Kroupa initial mass function instead 
of the Salpeter one wil l result in mass-to-light ratios a; 1 .7 times 
lower dMarastonll2005l) . 

7. Mass models 

We fitted the profile of the gravitating mass M{r) obtained by 
Eq. ([TJ with a mass model M mo a e \(r) that includes a luminous 
spheroid and a dark halo (Bertola et al. 1993; Kronawitter et 
al. 2000): 



M mode i(r) = M*(r) + M h (r). 



(6) 



In these objects the stellar surface brightness follows a Sersic 
profile with index n and effective radii r e given in Table [3] Once 
we deproject r e under the reasonable assumptions of spherical 
symmetry and constant stellar mass-to-light ratio, we obtain the 
spatial density p*(r) and mass profile M*(r) = J Am a pi,{r')dr' 
of the stellar component. 

Let us note that in NGC 7052 the X-ray mass measurement 
extends (only) to 16 kpc, corresponding to ^ 2r e . Not surpris- 
ingly, in this object the dynamical mass matches the spheroid 
mass over most of the probed region, and the contribution of the 
DM to the gravitating matter is small overall. Thus, to derive the 
structural parameters of the former component from the latter 
leads to results with sizable uncertainties. 

We assume as in spirals that the dark matt er is distributed in a 
spher ical halo with a Burkert density profile dSalucci & Burkeri 
|2000T) : 



P(r) = 



pari 



(r + r Q ) (r 2 + rl) 



(1) 



where the core radius ro and the central density po are f ree pa- 
rameters. The dark mass profile Mh(< r) takes the form: 



M b (r) = M { In (l + ^) - tan" 1 + i In 1 + Q 



(8) 



Ellipticals compl y with the inner baryon dominance 
parad igm like spirals dvan Albada et al.ll 1985t ISalucci & Persid 



The mass profile M(r) obtained from the above equations, 
with M mo dei (r, M*,po,ro) given by Eq. (O for our two ob- 
jects is fitted by minimizing the^ 2 . Estimated uncertainties are 
6 log M = 0.075 dex and 0.06 dex for NGC 7052 and NGC 7785, 
respectively. The results are shown in Fig. |5j the structural 
galaxy parameters, i.e., the spheroid mass, the mass, the core 
radius, and the central density of the halo are given in Table[3] 

The mass model with a Burkert halo and a Sersic spheroid 
gives a good fit to the distribution of the gravitating mass of our 
galaxies: almost all data are reproduced within the lcr error bars 



and^ 2 e 



1 . The transition region between the region dom- 



where M$ = 6.4por ( ,. 



inated by the stellar spheroid and the region dominated by the 
dark matter occurs at (2.0 + 0.3) r e . In both galaxies, we find 
that the stellar baryonic component d ominates the gravita tional 
pote ntial inside r e , in agree ment with lBertola et al.l ( 19931), and, 
e.g.. lFukazawa et ail d2006l) . At r e , the fraction of the gravitating 
mass contained in the stellar spheroid amounts to (76 + 5)% in 
NGC 7052, and (84+^)% in NGC 7785 where the lower spatial 
resolution makes the estimate more difficult. Then, the (M/L) 
profiles at small radii match those expected from the stellar com- 
ponent alone. 

We find that the stellar spheroid masses are M sp h = 1.1*2 x 
10 n M o in NGC 7785, and M sph = 3.6 + |^ x 1O U M in 
NGC 7052, where the quoted uncertainties are the fitting ones. 
The stellar mass fraction inside r e found here for both galaxies 
is quite consistent with earlier works (e.g., Gerhard et al. 2001; 
Thomas et al. 2009, and several recent SLAC papers starting 
with Koopmans et al. 2006). For NGC 7052 one by-product of 
the good agreement we find between the X-ray mass profile and 
the stellar one indicates that in this system there are no strong de- 
viations from the hydrostatic equilibrium condition, which may 
lead to an underestim ate of the cumulativ e mass, in agreement 
with previous studies dJohnson et al . 200§. 

In any case, we find that in NGC 7052 the X-ray derived 
mass profile is well fitted by the mass models. This is not ob- 
vious and in turn indicates that these measurements are very re- 
liable. Indeed, between 0.3r f and 0.8r f the mass profile of a de 
Vaucoulers spheroid has a unique distinct feature that we find in 
the data. 

Remarkably, these values agree very well with the above 
SED fitting estimates. Mass modeling leads in the B-band to a 
spheroid mass-to-light ratio of M sp h/LB = 6.3M /Lb o for NGC 
7052, and M sph /L B = 2.3M /L Bo for NGC 7785, with a fit- 
ting uncertainty of 30% (see Table 3). Note that the value for 
NGC 7785 is on the lower side for ellipticals, but it is confirmed 
by spectral fitting. In fact, in NGC 7785 the spectrophotometric 
estimate is M^/Lb = l-7M /L Bo (M s s ph = 8 x 1O 1O M ), while 
in NGC 7052 it is M^/L B = 6.9M /L Bo (M* ph = 4 x 1O U M ). 

Before we discuss the DM properties, notice that in 
NGC 7052 the dark matter dominates the mass distribution only 
at the outermost radii, for which we have X-ray data (1.5r e < 
r < 2.0r e ), which makes the estimate of its structural parameters 
uncertain. 

The dark matter halos show a cored density distribution. 
The core radii and the corresponding central densities take in 
NGC 7052 and NGC 7785 the values of 22.3/8.7 kpc and 
1.3 x 10~ 24 /1.0 x 10~ 23 g cirT 3 . However, these estimates are 
uncertain: even considering the mass profiles we infer from X- 
ray fluxes as error-free, the fitting uncertainties on these two 
parameters would be at the level of 50%. Furthermore, unlike 
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NGC 7052 NGC 7785 




r/r e r/r e 

Fig. 5. Left: Mass models of NGC 7052. The X-ray derived dynamical mass (filled circles with error bars) is modeled (thin black 
line) including a stellar spheroid (dotted line) and a Burkert halo (dashed line). Also shown is the best-fit solution for the total 
mass for models that assume a NFW profile to represent the dark halo (blue line) or the MOND framework (red solid lines). r e 
corresponds to 7.7 kpc. Right: Mass models of NGC 7785. r e corresponds to 5.6 kpc (see color figure online). 



for spirals where the dynamical mass is quite well determined 
from the rotation curve, here there is a non-negligible uncer- 
tainty in the (derived) mass profile that we fitted with the model 
6d\ogM{r)ld\ogr ^ 0.15, which contributes an additional un- 
certainty of ~ 30% on the derived values of ro and po. The com- 
parison of these values with those in spirals of the same mass 
must await more and more precise mass models; however, we 
can already claim that DM halo around ellipticals seem to be 
denser and with a smaller core radius, but the product poro is 
still, approximately, the same as the one found in smaller objects 
and for different Hubble types (Donato et al. 2009). 

Notice that in NGC 7785, similarly to what it is found in 
(some) other ellipticals, the X-ray derived mass profile in the 
very inner galaxy region, r < r e < 6 kpc, corresponding to the in- 
nermost 2-3 data points, cannot be reproduced by any reasonable 
mass profile, including a Sersic spheroid, a standard DM halo or 
a combination of both. The dynamical X-ray mass appears, for 
r — > 0, to be progressively smaller than the actual data. This 
discrepancy, probably caused by the contribution to the X-ray 
emission of i) non-thermal pressures or ii) multiple-temperature 
components (see Das et al. 2010) does not affect the main re- 
sults of the mass modeling, however. Indeed in the Burkert class 
spheroidal (BS) model the spheroid mass has been obtained by 
imposing that at 2r e , i.e., at a radius outside the troubled region, 
the spheroid and the dark halo masses equal the dynamical mass. 

Remarkably, we note that the resulting value for the spheroid 
mass is consistent with the photometric estimates. 

Finally, we checked the effects on the resulting best fitting 
values of the DM halo structural parameters (po and ro) of 1) 
a variation of the spheroid mass of +30%; 2) by considering 
the latter quantity as a free parameter; 3) including the two in- 
nermost data points. In all these checks, we found that the BS 
model fits the dynamical mass very well and the resulting val- 
ues of po and ro lie always within 30% those given in Table 3. 
The presence of a cored DM mass distribution relies on the fact 
that between r e and 3r e , i.e., well outside the troubled region, the 
dynamical mass shows an increase steeper than linear. 



We can estimate the virial masses and radii for NGC 7052 
and NGC 7785 by extrapolating the mass models: M v i r ^ 
8.6 x 10 12 M o , M vir = 5.2 x 10 12 M G , R vil = 476 kpc, and 
R v i r = 405 kpc. Indeed, in the present case the extrapolation 
of the outermost X-ray mass M(R ) to the virial radius is feasi- 
ble. The X-ray mass profiles indicate that the present data have 
reached, contrary to the "extended " available HI RCs in spirals, 
the region in which the circular velocity oc (M/R) 5 has a maxi- 
mum (Salucci et al. 2007). 

We are allowed to extrapolate M(r) because we know from 
weak-lensing measurements that DM halos around isolated 
galaxies are not truncated out to their virial radii (Mandelbaum 
et al. 2006) and that their outer profiles are consistent with a 
Burkert/NFW profile. The extrapolation uncertainties are at the 
level of 15% for NGC 7785, and 20% for NGC 7052. The un- 
certainty on M(R ) oc fjT(R ){3r~ l is also small, because the 
observational uncertainties on T(R ),/3, r c are at a level of 10% 
(Memola et al. 2009). 

These estimates of the virial mass from direct mass measure- 
ments are solid: these are two of the few cases in which this has 
been done so far. 

This confirms that big ellipticals live in halos two to four 
times more massive than those around big spirals. 

The baryon fractions, i.e., the ratios of stellar and virial mass, 
are 0.05 + 0.02 and 0.03 + 0.01, where the main uncertainty 
comes from how we extrapolate the dark mass profile out to the 
virial radius. The val ues found agree with the determinations of 
IShankar et al. (2006) obtained by correlating the baryonic mass 
function with the halo mass function. 

In summary, the fractional amount of the stellar component 
looks ( within a factor of tw o) similar to that found in the biggest 
spirals dSalucci et al. 2007]) and is in any case much smaller than 
the cosmological ratio of ~ 1/7. Given the large gravitational po- 
tential well of these galaxies, it remains to be understood where 
70%- 80% of the original cosmological baryons went. The dark 
matter mass profile in the central parts indicates a cored density 
distribution, although the estimate of the characterizing struc- 
tural parameters (central density and core radius) is uncertain. 
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7.1. NFW halo profile 

We further investigated the halo mass profile by fitting the dy- 
namical mass profiles with the NFW halo profile that emerged 
in the ACDM scenario (Navarro et al. 1997): 

m n m [log(l+^s)-(r/rs)/(l+r/r s )] 

M NFW (>) = M vk - — — , (9) 

log(l + c) - c/(l + c) 

with r s = R v i T /c. 

We tested this mass distribution, leaving c and M v ;,- as free 
parameters, which is necessary to obtain an acceptable value for 
the reduced;^ 2 ed . For NGC 7785 we leave the value of the stellar 
mass unchanged with respect to the one obtained above, because 
it is irrelevant for the total gravitational potential; for NGC 7052 
we leave it free. 

The model fits (see the dashed line in Fig. |5]l are very sat- 
isfactory for NGC 7052 (y 2 ed < 1) and less acceptable for 
NGC 7785 0r r 2 ed = 3.6). The resulting NFWS spheroid masses 
are reasonable and similar to the BS ones. In particular, in 
NGC 7052 the spheroid mass is 85% of the one obtained for 
the Burkert halo mass model. 

Moreover, in the latter case the resulting value for the con- 
centration parameter c — 33 is inconsistent with the ACDM pre- 
dictions of c — 8 - 9 for objects ~ 5 x 10 12 M o . This is not due 
to the mass modeling uncertainties, because the gravitating mass 
inside 30 kpc > 1O 12 M is much larger than the spheroid mass, 
therefore coincides with the halo mass, which implies that this 
object displays an intrinsically cored mass distribution. 

As a counter proof, the standard NFW DM halo of \0 l2 M^ 
has c - 9 (e.g.. jBertola et al.1 1 1993b IKronawitter et ai]l2000t 
iKlvpin et al.ll2010t lTortora e t alJl2010h and it is completely in- 
consistent with the data. 

Note that halo models with c 3, which is equally incon- 
sistent with ACDM p rediction, but in line with the results of 
iFukazawa et al.l d2006l) . have a ^ 2 ed similar to cored best-fit val- 
ues. 

Changing the values of the spheroid masses or considering 
the adiabatic contraction process for both galaxies for the model 
with the NFW halo and Sersic spheroid (NFWS) gives a worse 
fit. The NFWS model vs observation discrepancy is very dif- 
ferent with respect to the BS one. Firstly, it extends outside the 
(possibly) troubled region, out to r < 2.5r e < 15 kpc. Secondly, 
the discrepancy is also found in a region where only one com- 
ponent is present, the dark one, which completely dominates the 
potential with a known density profile. Furthermore, the (very 
high) value of c does not depend on the innermost values of the 
dynamical mass, but on the fact that from r ~ r c ^ 1.5 r e ~ 8 kpc 
out to the last measured point at ~ 45 kpc the dynamical mass 
(completely dominated by the dark component) increases per- 
fectly linearly with the radius. Independently on how we treat 
the first two data at r < r e , the NFWS mass model for NGC 
7785 must have c ^ 35. 

7.2. MOND 

We realize that our objects are in the Mondian regime for r > 
5 kpc. Our two objects are then well suited for testing the MOND 
paradigm, in particular the case in which it is seen as a modifica- 
tion of the inertia law, able to account for the "missing light" 
in galaxies. In this framework the relation between the den- 
sity of the gravitating matter and its relative acceleration, when 
the latter is < «o = 1.3 x 10~ 8 cm s~ 2 , is very different from 
the Galilean-Newtonian case. As a consequence, for a baryonic 



mass distribution Mb(r), we measure (e.g. by applying the hydro- 
static equilibrium condition) a different/larger dynamical mass 
MmondW: 



MmqndOO = ^M h (r) 1 1 + 



1 + 



4.1 x 10 9 M o 



lkpc 



0.5 i 



■(10) 



Therefore, we test MOND by fitting M(r) given by Eq. (O 
with Mmond(/); leaving the mass of the spheroid as a free pa- 
rameter. We can anticipate that MOND will have difficulties ac- 
counting for the DM of early-type systems. They have just one 
baryonic component, the Sersic spheroid, whose mass converges 
well inside the region for which we have the dynamical mass. 
This creates a mass discrepancy that steadily increases with ra- 
dius as a power-law. In spirals MOND must account for a very 
similar behavior of the dynamical mass but it has two additional 
advantages: 1) there are two baryonic components - the stellar 
and the HI disk, and 2) the mass of the latter does not converge 
inside the region mapped by dynamical data. 

MOND shows a tension in NGC 7052, and fails the test for 
NGC 7785 (see Fig. 0). The value x 2 Kd - 2 of this mass model 
is fairly high. However, even if we consider the MOND mass 
model satisfactory, the resulting best-fit values for the spheroid 
mass-to-light ratios (M/Lb - 7 and 15) are inconsistent with the 
colors of these galaxies. 

8. Conclusions 

We study the mass distribution of two isolated ellipticals by 
means of their X-ray emission and stellar photometry. We find, 
and this is one of the few cases at present for early-type galax- 
ies in such suitable environment, that a dark component must 
be present for r > r e , because r e is the half-light radius of 
the stellar spheroid. More specifically, we find that inside r e 
the stellar component accounts for $ 70 % of the gravitat- 
ing mass. This result is not new (see, e.g., IKronawitter et al.l 
l2000tlNagino & Matsushital2009L where 22 objects with XMM- 
Newton and Chandra X-ray luminosities are mass-modeled). 
However, our finding is at o dds with the claim of s ubstantial 
proportions of DM inside r e dHumphrev et al.l 12006) and with 
the idea that the well known tilt of the fundamental plane of el- 
lipticals is due to systematic variations of the DM fraction inside 
r e with galaxy properties. Indeed, our (and other) mass models 
of ellipticals imply that the central stellar line-of-sight velocity 
dispersion has very little contribution from the dark component 
and essentially weighs only the stellar spheroid, whose struc- 
tural properties (mass-to-light ratio and half light radius) there- 
fore must be considered responsible for the observed tilt. 

In our two objects the total mass does not increase 
with radius as a (uniq ue) power-law, as suggested by 
iHumphrev & Buotel (120101) : dlogM/dlogr decreases from ~ 2 
in the innermost regions mapped by X-ray emission (~ 0.3r e ) to 
~ 1 in the outermost regions (~ 2r e ) , a behavior often found 
in the mass structure of ellipticals (IKronawitter et al.l 120001: 
IFukazawa etaD 120061) . The mass profiles of NGC 7052 and 
NGC 7785 are easily and naturally reproduced by a combina- 
tion of a Sersic spheroid and a Burkert halo, where the former 
is much denser and the latter is much more massive. The fits 
are good (unfortunately) for a large range of values of the struc- 
tural parameters, i.e., the central DM density, the core radius, 
and the spheroid mass. This agrees with previous results: ellipti- 
cals show a variety of mass profiles, each of them reproducible 
by a quite large variety of different mass models. Thus, in our 
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Table 3. Structural properties. Columns: (1) galaxy name; (2) last measured point, in arcsec; (3) half-light radius in kpc; (4) Sersic 
index; (5) mass inside the last measured point, in units of 1O U M ; (6) spheroid mass in units of 1O U M ; (7) mass-to-light ratio of 
the spheroid in B-band, in units of M /Lb o ; (8) dark matter Burkert halo core radius in kpc; (9) dark matter Burkert halo central 
density in units of l(T 24 g cirT 3 ; (10) Burkert halo virial mass in units of 1O I2 M ; (11) concentration parameter c of NFW halo; 
(12) NFW halo virial mass in units of 1O 12 M ; (13) Burkert halo virial radius in kpc; (14) baryonic fraction. 



Source 


lpt 


r e 


n 


Mi pt 


M sph 


M sph /L B 


ro 


Po 


M vir 


c 


M VNFW 




./baryon 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


NGC 7052 


16.1 


7.7 


4 


4.0 


3.6 


6.3 


22.3 


1.3 


7.4 


10 


8.6 


476 


0.05 


NGC 7785 


32.1 


5.6 


3.8 


9.0 


1.1 


2.3 


8.7 


10.2 


3.6 


33 


5.2 


405 


0.03 



objects we find no evidence for the claim that the mass distri- 
bution in ellipticals is subject to a "cosmic conspiracy", i.e., that 
all of them would show the same (unexplained) profile, implying 
just a relationship rather than the claimed fine-tuning in the DM 
and luminous matter structural parameters. 

As in spirals, NFWS mass models have some difficulties re- 
producing the data: the best-fitting values of the x 2 , the mass 
to light ratio of the luminous component, the halo concentration 
parameter are out of the expected range. However, in ellipticals, 
a NFW profile is not necessarily the correct present-day profile 
for a ACDM halo as well. The baryons largely dominate the in- 
ner potential and thei r infall/contraction may have modified the 
DM distribution (e.g jKlvpin et al.ll2002l) . 

Finally, our mass deter minations are in tension with the bary- 
onic Tully-Fisher relation (iMcGaugh et al.ll2000h . according to 
which in cosmic structures from clusters to dwarf spheroids 
(dSphs), the circular velocity V = GM(r)/r 5 and the bary- 
onic mass Mb scale as log Mb = 4 log V + 1.61. In our objects 
we obtain a reliable estimate of the stellar mass that coincides 
with the baryonic one. Moreover, for our objects the mass pro- 
file does increase linearly with radius in the outer regions, so 
that the "circular velocity" is radially constant with radius and 
can be easily extrapolated out to its virial value. We find that 
for NGC 7052 (V, M b ) =( 330 km s -1 , 3.6 x 10 n M o ), while for 
NGC 7785 (V, M b ) = (380 km s" 1 , 1.1 x 10 n M o ). These objects 
have baryonic masses that fall short of the predicted value from 
this relation by a factor ~2 and ~8. 

By investigating two isolated ellipticals by means of their 
X-ray and optical emissions we find that 1) they have a consid- 
erable amount of DM, but only at radii > r e ; 2) their DM halos 
have a shallower distribution than that emerging in the ACDM 
scenario; 3) these objects are in the Mondian regime, but do not 
agree with its predictions; 4) these objects do not comply with 
the baryonic Tully Fisher. 
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